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'KATIONAL ADVISOR! COMMITTEE FOR AEROWAUTICS 


TECHNICAL NOTE NO. 1170 

ON SUBSONIC COMPRESSIBLE FLOWS BY A METHOD OF CORRESPONDENCE 
I - METHODS FOR' OBTAINING SUBSONIC CIRCULATORy COMPRESSIBLE FLOWS 
ABOUT TWO-DIMENSIONAL BODIES 
By Abe Gelbart 

SUMMARY 


By means of the general eolutiraia of the hodograph equations for 
coQpresaible fluids, certain solutions corresponding to solutions of the 
hodograph. equations of an incompressible fluid are used to find flow 
patterns of conqjresslble fluids. When the adiabatic equation of state 
is used, only a general method is outlined. 

The method appears to lead to the solution of the problem of sub- 
sonic flows with circulation around arbitrary bodies, as the methc^ of 
Theodorsen does for Incompressible fluids. A second paper, part II, 
illustrates the method for some given bodies. For the linearized equa- 
tion of state, the results obtained Include some of the results of Von 
Harman and Tslen, as well as some of the recent results of Bers. The 
method can be used for flows with circulation as. well as without circu- 
lation. 


INTRODUCTION 


It is well known that the nonlinear con 5 )re 8 slble— flow equations in 
the physical plane can be reduced to linear equations which have the 
form of generalized Cauchy-Riemann eq\iations by the change of the inde- 
pendent variables in the physical plane to the independent variables of 
the hodograph plane. The successful technique applied to the solution 
of the incoirq)resslble— flow problems around given bodies was achieved 
only because the theory of analytic functions had been developed pre- 
viously. It would, therefore, be suspected that similar results for the 
theory of con^jressible fluids could be obtained if there were developed 
a theory corresponding to the theory of analytic functions for the gen- 
eralized Cauchy-Riemann equations representing the flow of a compressible 
fluid in the hodograph plane. 
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Develqpment of such a theory has "been, partly achieved by the theory 
of E-inoiiogenlc functions (references 1 a n d 2} . This theory has not 
yet attained the perfection that the theory of analytic functions has. 
However > the S-inonogenio Taylor series and the s-monogenlc, Laplace 
transforms JJiay be partly able to overcome the weaiaaesses of this method. 
Some eiTorts in this direction sire being made in this paper. . . 


Even with the complete theory of S-monogenic functions one of the 
chief drawbacks of the approach lies in the difficulty of transforming 
the solutions from the hodograph plaiw to the physical plane. 

The method of the theory of correspondence has been very briefly 
outlined in reference 3* The procedxire is to obtain the flow of an in— 
can 5 )reeslble fluid around some given closed body, to transform the com- 
plex potential of the flow to the hodograph, and then to obtain the 
particular S)-monogenlc function of the infinite set of solutions which 
possesses the desired prqperties. 

It may seem, at first, that this technique is too general to be of 
practical use, but there already exist some interesting results of the 
application of this method. Fiirthermore, this technique is certainly 
not new. It has been used with much success by Chaplygin, Ton Kar man , 
Tslen, Bergman, Here, and others. (See references 4 to 12.) However, 
much is believed to be new in the specific use of this method when 
eriploying 2-Jncnogenic functions. When the method is applied directly 
to that of a source and a sink of an incong>ressible fluid, the method 
yields .a source and a sink of a compreseible fluid. 

It is the alia of this paper to elaborate on the correepondence 
method in a general way and to apply it in particular to the flow of a 
coittpres Bible fluid under the linearized equation of state. It should 
be mentioned that this method has already been applied, under the as- 
sumption of the linearized equation of state, by others, notably by 
Chaplygin, von Kerman, Tslen, and Bers. 

By dealing with the method in all its generel.ity, a formula is ob- 
tained involving an arbitrary analytic function. By choosing particular 
values of this arbitrary function, of which the derivative is regular in 
the exterior of a region including the origin, the particular foimulas 
of Tslen (reference 6) eind Bers (reference 8) are obtained. Once the 
arbitrary analytic function Is determined for a given flow, velocities, 
stagnation' points, and so forth, can be readily computed. In a second 
paper by Bar'fcnoff and Gelbart (reference 13 ), some flows around given 
bodies are calculated. 

This Investigation was carried out at Syracuse Iftiiversity xmder 
the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeromutics. This report was submitted in July 
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SYMBOLS 

k, a, p constants 

p pressure 

p density 

q. velocity 

Pi# <3.1 particular vaJLues of the (juantlties 

o subscript referring to the state of the fluid at rest 

u horizontal velocity 

V vertical velocity 

6 angle that velocity vector makes with horizontal 

Q velocity of undisturbed stream 

*00 

q distorted velocity 

a velocity of sound 

w conplex variable in distorted hodograph plane 

a complex variable in hodograph plane 

z complex variable in physical plane 

M Mach number (q/a) 

qp, $ potential functions 

\j; stream functions 

Q complex potential (<p + l\i/) 

^ complex variable in auxiliary plane 

ratio of specific heats 
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K constant depending on velocity of undisturted stream 

K.( 

Gr(^) complex potential of an Incong^ressible fluid In auxiliary plane 

r circulation of an incompreesi'ble fluid 

T an analytic function of the velocity q. 

(n) ~(n) 

W ,iW generalized complex powers 
generalized real powers 

\ 

generalized differentiation 
• generalized multiplication 

13100111181. coefficients 

E generalized complex eacponential function 

c,c* generalized real cosine function 

B,s* generalized real sine function 

L-j, generalized Laplace transform 

( ) coii 5 )lex conjugate of ( ) 

absolute value symbol 

COMPHESSIBLE ELOVTS UH12EE THE ASSUMPTION OE 
THE ADIABATIC EQUATION OE STATE 

The four basic relations for the potential flow of a steady two- 
dimensional fluid that will be assumed are: 

Equation of state, p = (l) 



where 
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P 

P 

7 

k 


presBiire 

density 

ratio of the specific heats 
constant 

Bernoulli *s eq.uation. 



dp 


= constant 


P 


( 2 ) 


The continuity equation. 


div 



and the circulation eq.uation (for irrotational fluids), 

curl = 0 . 

The Telocity of sound is given hy 


(3) 


W 


dp 


(5) 


From. eq.uation (5) the Bernoulli eq.uation can he written in the form 


i! 

2 



= constant 


(6) 


which often is more convenient. From these assiuigitions on the fluid, 
the first two equations give rise to the relations fPr density, pressure, 
the velocity of sound, and so forth, in terms of the velocity only. The 
third and fotirth equations give rise to the equations of motion of the 
flow. 


The subscript zero on the variables’ p, p, and a will indicate 
the particular value of the variable at a stagnation point . It will be 
assumed throughout that 
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(7) 


This la equivalent to introducing the dimensionless variables p/Pq 
and q/ag. 

From equations (l) and (2), with the aid of equation (5)^ the fol- 
lowing well-imown relations are established; 

sfi « - [ (7 - l)/2 ] q2 (8) 


P 


= p 


o 


P =* P 


0 


S3 


Pq 




_af 

2 ao^ 


) 


_2_ 

7-1 


+ 


_ 7 

y - 1 if N 7-1 
2 a^y 


1 

_ y - 1 

p o 2 y 

+ 7.. - 1 a!V ^ 

2 a^y 


M® 


q £P 

P dq 



(9) 


( 10 ) 


( 11 ) 


These quantities are all given in terms of the single v«iriable q. 


From the basic relations (3) and (4), and the fact that the flow is 
potential, there exist two functions <p, the potential function, and \1/, 
the stream function, that satisfy the equations 


'Px 



9, 





( 12 ) 


J 
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where the subscripts indicate partial differentiation with respect to the 
variables indicated. The independent variables 2 and y are the coor- 
dinates in the physical plane of the position of the particle of the 
fluid. System (12) represents the flow of a compressible fluid. 


Let u and v be the horizontal and the vertical velocities, re- 
spectively, of a particle of the fluid at a point (x,y) . Then 



( 13 ) 


where q. is the magnitude of the velocity q., and 9 the angle that q. 
makes with the x-azis. 


From the physical definitions of the potential function cp and the 
stream function \J/ , it follows that 


dcp = u dx + V dy 


( 14 ) 


d\J/ = — pv (3x + pu dy 


(15) 


Since dcp 5 - cp^ dx + cpy dy and d\j/ = \]/^ dx + \j/y dy, system (12) can bo 
obtained directly by ccmgjaring these eq.uations with eq.uations (l4) and 

(15) . 

From equations (14) and ( 15 ) it follows that 

dcp+ i — d\!/=udx + vdy + i(— v dx + u dy) 

P 


= (u — iv) (dx + i dy) 

= qe~^® dz 

where z = x + iy. Thus it follows that 

19 ^ 

dz = 2 — ( dj» + ^ d\j/ 
q N P 


( 16 ) 


In the hodograph, the coordinates are the polar ^ coordinates (q,0). 
■Equation ( 16 ) is in the nattire of a transformation from the hodograph 
coordinates of the flow to the physical coordinates (x,y) . It is this 
fact that makes equation (I 6 ) of fundamental importance for the approach 
used in this paper. 
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I^on eliminating first one and then the other of the two iinknown 
functions, cp and \(/, from equations (12), the second-order partial 
differential equations 

(<^X + W y ' “ 

(5'l'*)x- (I*)')?'® 

aie obtained 

I'i'om. equation (l4) It is deduced that 

q = + <Py^ (^9) 


so that p is a function of the unknown cp. Equations (l?) and (l8) 
are therefore nonlinear. It is precisely this condition that has led 
previous investigators to ti’ansform the flow equations from the independ- 
ent variables (x,y) in the physical plane to the Independent variables 
(q,6) in the more geometrically conqjlicated holograph plane. As will be 
shown, these geomotrio conqilications can be oircxuaVonted after the lin- 
earized equations in the holograph variables are used to advantage. 


The flow equations in the holograph variables have been derived 
from equations (3.2) by many previous writers. For the eeike of ccmplete— 
ness, however, an outline of a derivation will be presented here. (See 
reference 10.) 


Equation (l6) is first differentiated with respect to 0: 


hd q V50 p hd J 


and then with respect to q: 


(20) 


^ 1 + i (21) 

3q q '^q p 8q' 

Equation (20) is differentiated with respect to q; 


S^z i0 

r ^ H 

1 

< 1 1 

1 , 1 1 i ^ 

w ®* © 

Sq30 

q2 30 

<iq 

^pq heJj 

[ q Vdq 30 p 8q 30/^ 


( 22 ) 
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and eq.viatlon (2l) is diffearentiated with respect to 0; 


Jfz_ = i + i 1 ^ 1 S-- + 

30 Sq. q. \ 3q p 3q. / q. \ 30 3q. 


1 i 

p 30 3q. / 


(S3) 


Since the leftHiand sides of eq,uations (22) and (23) are eq.ual, it fol- 
lows that 


a® Ug^"p3a; ® qt^ 30 " da Uq. J a© 


(24) 


The flow eauations in the hodograph variables are obtained by equating 
the real and the imaginary parts of equation (24): 


^ = a M 

30 p 3q 

> 

3q dq \ pq/ 30 


(25) 


From equation (ll) 


§P = _ £ 

dq q 


It follows that 


i_ ^ L. 

dq\pq/ p^q dq pq^ 

= A- 14^ = -A. (M^ 1 ) 

pqL^ pq.^ 

Equations (25) can now be written in a more suitable form: 


3q) _ a_ 

30 p 3q 

^ = _ 1 - ^ 

3q pq 30 


(26) 
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The fact that M and p • are functione of g only makes eguatlons 
(26) linear. These eguations can he handled mathematically with much 
greater ease. than can the nonlinear eguations (12) in the physical plane. 
The theory of E-fflonogenic functions has been set up for just this pur- 
pose. Thoiigh'the solutions of eguatione (26) that are obtained give the 
complex potential of the flow in terms of the velocity of a particle of 
the fluid rather than in terms of- the position of- the particle, the fun- 
damental transformation, eguation (l6), enables these solutions to be 
transformed back to the physical plane. The fact that the right-hand 
side of eguatlon ,(l6) is an exact differential facilitates greatly the 
process of carrying out this transformation. 


SOLUTIONS OF THE COMERESSIBLE FLOW EQUATIONS AS FUNCTIONS 
OF THE EODOGRAPH VARIABLES 


Equations (26) can be written in a more genoreil form for the purpose 
of mathematical treatment: 


dp _ 
be ~ 

dq 


Ta(q) 
- Ta(q) 


d\(/ 

dq 

de 


(27) 


where and T are positive analytic functions of q. The condition 

that Tg be positive is equivalent to assmning that the flow romalnrf sub- 
sonic. Many of the mathematical results that will be preoented here hold 
also when Tg becomes negative, so that flows involving supersonic ve- 
locities can be treated. This has been carried out elsewhere for very 
special flows. (See reference 3*) 


Consider the function represented by the line integral 


F(c) = 0(0, q) + 1'$ (0,q) 

pd0 -T^ (q)\i/ dqj + if jVd0 + - — ^ 9 dq^ ( 28 ) 


0;q. 

r*» 




'o^*lo 


Of^O 


where Qq, q© are fixed values of 9 and q and a = 0 + iq. By 
virtue of equations ( 27 ) a sinple computation verifies that each of the 
two integrands in the right-hand side of equation ( 28 ) is an exact differ- 
ential . 
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This establishes that F Is a function of the point ( 6 >q.) and 
has the same value when the Integration is taken along any path, p 2 rovided 
prefer regard is given to singularities. 

The second important fact about F( 6 ,q.) is that its real part $ 
and its imaginary peirt $ again satisfy the system of equations { 2 '( ) . 
This affords a method of generating particular solutions of system (27 ) . 


It is trivial that when <p = 1 and \1; = 0, this pair is a solu- 
tion of eq.uatlons (27) . When this solution is substituted into eq.ua— 
tlon ( 28 ) 


$ + i^‘=* ¥(ct) 

=» S + i 



(29) 


If the real and the imaginary parts of eqxiation (29) are again substi- 
tuted into eq.uation ( 28 ), the next solution gives 


W 


( 2 ) 



(30) 


The superscript 2 in parentheses indicates that the function has the 
nature of a power, as is observed when ”’’1 *’’’2 = 1* Here equation ( 30 ) 
reduces to (9 + iq)^ = a^. 

The notation F(aQja), Cq » 0q 
are the lower limits of integration 

= (9 - + 2i (9 - 


and corresponds to the function (a 
usually written F(a).' 


+ iqQ, indicates that 9 q and qg 
in equation ( 28 ) . For example. 



— Oq) • When = 0, F(cTQja) is 


Let 



1 


and 
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n odd 

(31) 


n even 


If this operation is repeated n times the resulting solution is 

= (0 + iQ)^^^ 

n 

u = o 


n times 


(t) = 



where C„ ,, is the vth "binomial coefficient. 

11^17 

If another trivial solution is taJsen as a first solution, say, 
cp a 0, \(/ a 1, and substituted in equation (28), it follows that 


F(a) a i.W(a) 
= iW(o) 


1(0 + i/ Tg dq) 



The symbolic notation i.W is used to represent a genoradlzatlon of the 
concept of- multiplication. It will be seen that tMs notation is rather 
useful. When this process of substitution is repeated 



and 
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n times 



Again, if this operation is repeated n times the resulting solution is 


= 1 


n 

I 

u=o 


■'n,u 




(36) 


Finally, if the trivial solution cp = a, \]/ = P is tahen, a and p 
heing constants, nr-repetltlons of the operation (28) yield tiie solution 



= (a) + p 



( 37 ) 


where a = a + ip. 

Coicplex— valued functions of which the real and imaginary parts sat- 
isfy the system ^ 26 ) are said to be 2-monogenlc functions. Thus, 

and = IW^^^ are s^Jionogenlc functions. These are more 

specifically referred to as fomnal powers, since when T ^ * Tg = 1, they 

reduce to = (6 + iq.)^ and lo^ = 1{S + lq.)^« The formal product of 
the two solutions a and is a*W^^^, which Itself is a solution. 

It should be noted that the ordinary product aW^^^ is not a solution 
of the system (26) when a is coniplex. 

As S -multiplication of a constant end a formal power was intro- 
duced, so may E— integration be introduced. The S-integration of the 
2-monogenic function? f (a) = <p + i\Jf, is defined as 

F(o) =0 + 1$ 



f (cr)d^ a 


(38) 
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The ^ refers to the coefficient aiatrix of system (27) . 


) « 


1 T. 


U T, 


(39) 


And E-differentlatlon of the 2-monogenic function, f(c) » q> + ij|/, 
is defined as 

V, , d f(o) 
f(a) « — ^ 


~ 4bc + i\{/x 

From the relations of equations (27) 


('to) 


dr ^ 


= '’'i'lty-f-j«Py 


(41) 


Again, when Tj. *» Tg = 1^ E— integration and 2-diffs3rentiatlon 
Induce to ordinary coaplex integration and differentiation. 

A more elaborate definition of E-^integration and differentiation 
has been given in reference 2. It has been shown in reference 2 that 


and 


dg F(o) 

dj,a 


f(o) 


(42) 


y^F wdj. 


P(ct) 


(43) 


Thus, ^“integration and S-differentiatlon are inverse processes . Also, 
the S-derivatlve of a E -monogenic function is a S-monogenio function. 


From the preceding definitions it can be verified by direct coacpu— 
tatlon that the E— integration and the E -differentiation of the formal 
powers follow similar rules to those of ordinary integration and differ- 
entation of the powers of a con^plex variable; more precisely. 



a*W 


(n+i) 


n + 1 


( 44 ) 
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and 






na*W 


(n-i) 


( 1 ^ 5 ) 


Since system (2?) Is linear, the sum of two solutions is a solution. 
Hence, a formal polynomial of the nth degree. 


f (c) = ao + aj_*W + . . » + aji*W^ ^ , aj^ ^ 0 

is a 2 -monogenic function. The formal power series 

00 

f(a) = ^ (47) 

n=o 


represents a E -monogenic function, p2X>vided the series converges uni- 
formly and absolutely in some neighborhood of Co- 


lt has been shown (reference 2) that every complex solution of 
system (27) can be represented uniq.uely by a formal power series of the 
form of eq.uation (47) - (This statement is true only when Tj_ Tg > 0, 
that is, in the elliptic case or subsonic flow.) The coefficients are 
given by the formula 

n 


Sn 


n! 


d£f(c) 

d^CT 


cr = Cr 


(48) 


Since every solution of system (27) can be represented, at least 
in some neighborhood, as a formal power series and since every coBg)lex 
potential in the holograph variables of a compressible fluid is a solu- 
tion of system (27), it would appear that it remains only to choose that 
particular solution which is the desired flow around a given body with 
prescribed conditions . This approach, however, still has some very 
serious difficulties. 


The formal powers presented previotisly are solutions of system (27) 
in closed form and are valid throughout their regions of regularity. 

The first powers W(a) and i-W(a) are laiown to represent a compress— 
Ible vortex, respectively. The higher powers appear not to represent 
any flows of interest. Linear combinations of the powers, that is, 
formal polynomials, have not yet been studied sufficiently to determine 
whether they represent flows of interest. The solutions represented by 
formal power series are first, not in closed form and, second, not valid 



l6 


If AD A TN No, 1170 


throughout their regions of regularity. Furthermore, since in aerodynaat- 
ics the majority of solutions of Interest are those that are regJilar 
throughout the exterior of a closed region, solutions corresponding to 
the analytic functions o~^ must he obtained in either closed form or 
in a series that can he extended throughout its region of regularity. 

What is most desired is a theorem analogous to that of Laurent in eooaT- 
lytlc functions, for 2-®ionogenic functions. Such a theorem, if it can 
he found, might quickly lead to the solutions of the major prohlems in 
compressibility , 

It is hardly hoped that the coa^ilicated problem of compressibility 
can be solved in any einqjle way. Any results in this direction, there- 
fore, are of Interest. 

Some progress along these lines has been made by Bers and Gelbart' 
(reference 2) by extending to E-Jnonogenic functions some of the results 
on the Laplace transform. These give, in closed form, solut^ns of 
system (27) that are different from the foiinal powers &>\fw , and in 
half planes correspond to the analytic fvmctions cT^. 

Consider the function 

00 

a-E(Coja,a) = ^ ^ _ (49) 

n=o 

where as before a « a + 10 (see references 1 and 2) and a is real. 

For the sake of brevity 


l.E(Co;a,cy) = 

E(Oja,cr) = E(a,c) 

This function is termed the ”2r-esponentlal function" for obvious 
reasons . ’ . . 

A simple confutation shows that 

E(a,a) = e®^[c(a,q) + is(a,q)] (51) 

and 

l*E(a,a) = ie®® [c*(a,q) + l8'^(a,q)] {^2) 



where 
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n=o 


V (-1)“ „.(==) 




n=o 


s(a,q.) 


s-^ (a, q) 


y iril^ 


^ai+1 


Q 


(an-ti) 


n=o 

03 


(2n + 1) ! 

.a an+i 


y (ziT^,. 

(en + 1)» 


(a) 

(an+i) 


(a) 


n=o 


(53a) 


(53b) 


From til© definition of E(a,cr) and that of S -differentiation, it 
is clear that 


^ 7 : 

d y O' 


[E(a,c)] = a E(a,a) 


Further properties of functions c, s, c* , and s* are 


m 


s(a,0) = s’^(a,0) = 0 
c(a,0) = c*(a,0) = 1 


and 


s* (a,<i) 
s*‘ (a,q.) 
c ‘ (a, q.) 

c*» (a,q) 


g 

Tx(<l) 

a T 2 (q.)c--'(a,y) 

-a ''2(q.)s(a',qL) 


= -a 


8* (a,q.) 
■^i(a) 


J 


(55) 


(56) 
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tti0 differentiation teing with respect to q,. These functions may be 
regarded as generalizations of the trigoncanetric functions and have many 
propearties aoAlogous to those of the trigonometric functions . For 
example, they^Satisfy Sturnt-Liouville differential equations and possess 
such properties as 


o(a,q)c^(a,q) + e(a,q)s*"(a,q) h 1 


( 57 ) 


How consider the function 


f(d) 


= [F(a)] 


CD 



0 


E(-a,o)da 


(58) 


where F(a) = Fx(a) + iF2(a) is a coa5>lex-valued function of the real 
variable a, Fi(a) and Fa(a) being real functions. This is called 
the S -Laplace transform of F(a), for when = Tg = 1 the Integral 
in equation (58) reduces to the ordinary Laplace transform of the func- 
tion F (a) . 


Equation (58) may be rewritten in the znore convenient form 


f (a) = q>( 0 jq.) + 


Fi (tt)E cf )cLtt + j ^2(^0) 0 t» cj )3 dot 


0_03 


» / e~®^ c(a,q)Fi(a)da ~ i / e^"^ s(a, q)Fi(a)da 


r-ClB 


/ 


+ i / o~^ c* (a,q)Fg(a)da + / e"*^ s*(a, q)Fg(a)da 


-O 0 


Then 


( 59 ) 


00 00 

0 (a,q) = J* c(a,q)Fx(a)da + s* (a,q)Fg(a)d<x. (60) 


and 


o 0 

00 00 

OC0 _\ t !, CC^ 


\ 1 / ( 0 ,q) - — I s(a,q)Fx(ct.)cLa + / e"^ c'*' (a,q)Fg(a)da (6I) 


Two facts may be demonstrated about the function in equation (58) : 
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1. If 3 is the ahsclssa of convergence of the Laplace Integral 

F(a)da (62) 

then the generalized Laplace integral 

F(a) • E(-a,a)da (63) 

converges uniformly in the semi— infinite strip 




0 > 0q>3 

< q. < 42 


2. Integral (63) represents a S-fflonogenic function in the open 
half plane of convergence of integral (62) . 

Because of theorem 2 , expression (63) represents a new class .of 
solutions of system (27) in closed form. Hence, for every fiuiction 
F(a) for which integral (63) converges (from theorem 1 this class is 
known to he wide), expression (63) is the complex potential of a com- 
pressihle flow in the holograph coordinates. 


For certain functions F(a) it is easy to show that the formal 
series expansion of function 


f (a) 


GO 


= / F (a) •$(-<[,, a) da 


converges in a smaller region than does function f(a). The E— Laplace 
representation of a function may thiis he regarded as a means of contin- 
uing analytically, a function that is given in the form of a formal power 
series . 


Because 


_1 

a 


-oc , 
e da 


(64) 


in the right half plane, it seems natural to define the negative S — 
power function hy ^ 

» J" S(-a,a)da (65) 
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Since E -differentiation under the integral sign is permitted, rep- 
resentations for the higher negative S -powers can he obtained. By E— 
differentiating each side of equation (65) n times and using equation 
(5^), it follows that 





n 



da 


( 66 ) 


In view of equation (66) it seems convenient to define the negative 
E— powered functions by 


W 


(-n) 


n 


CO 


t E(-a,o) 


rt-i 


a 




(n - 1) I 


da 


(67) 


Generalizations of other special functions of a complex, variable may be 
obtained in a similar way. 

In order that functions as defined by equation (67), pday 

a similar role to that of the Inverse powers it would be necessary to 
show that they possess a pole of the nth order at the origin and be . 
regirlar everywhere else. This has not yet been shown. In fact, wV 
is defined only for the right half plane, Ela > 0. It can be defined 
in the loft half plane, Elcr < 0, by 

^( n) _ r 2(a,o) — — da (68) 

J (n - 1) I 

o 


On the imaginary axis other than at the origin it can be defined as the 
limit of the function as a approaches a point on the imaginary axis. 

It has not yet been proved that defined in the right and the 

left half plane by equations (67) and (68), respectively, is continuous 
across the imaginary axis. If its truth is assumed, some progress might 
be made in the study of compressible flows. 

In order to study the flows that arise from these solutions, it is 
necessary to transform them from the hodograph back to the physical 
plane. 


Consider the solution 

cp + 1\1/ ® 


i^ C 


n.r 


n=o 


0n-r ^(r) 


(69) 
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The real part is 


t 


<P 


n-ar (2r) 


X .r n-ar 

(— 1 ) ^XifST ® ^ 


(TO) 


r=o 


and the Imaginary part is 


[v] 


n— ar-i (arf-i ) 
Q 


( 71 ) 


r=o 


where the "brackets in the upper limit of the summation indicate that the 
integer chosen is the smaller one nearer to the number within the 
brackets . 

After substituting (jO) and (71) into equation 

e^® , i 
dz = (dcp + - d\j;) 

q. P 


,10 r 1 1 

<pg d 0 + 9 da + - (\j /0 dS + tg da) 

qi L ■ p ■ • . J 


(72) 


and integrating, it follows that 

an-s+>(-i)^+i' 


z = e 


{i I [(-l)*"(n-2r)0„„, 

' 16-1 („-ar -l)! 

r* / Mi+n * — ' - J3. •“ "1 

1 an- 3 +(— ] s=o 

(-!)"'(» - Sr - l)C^,srn 


L 


I 


r=o 


-] 


-ar— 1 — s 
s 


n-ar-a 


y / \ r „n— ar-s-*a -1 

iS Q 1 ( 73 ) 

(n - 2r - s - 2) ! _, j ^ 


B=0 
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Eq^uation (73) together with equation ( 69 ) cem he regarded eis the 
parametric equations of the oompresalhle flow of cp + 1\|/ In the phyal— 
ceil plane. 

Similarly, for the solution 
cp + 1\1/ = 1W^“^ 

= 1 T 

r=o 


en-3+( 


= <1 


z ~ e 


4 


^ I" (-l)'=’(n-2r)Cj,^2r . 


r=o 


'gn-s+(-i)^+^ 


+ 

PI 


y 


r^o 


n-ar-i 


— - , ^ , s f _n-23>-e-i -i 

Y* ^ 6-1 (n - gr - 1) I e 

— (n - 2r - B — 1 ) 5 


B=0 


■(-if (n - 2r - 


n— 2 r ~2 


y iS (n - gr - 2)1 e 


li—sy—B—z 


s=o 


(n - 2r — s — 2 ) I 




(75) 


By taking linear comhlnatlons of the formal power solutions, other 
flows can he obtained In the physical plane In parametric form. 


If the generalized Laplace Integral solutions are taken, then 
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q> + i\l/ = / F(a) •E(-a,a)d0 


and 


00 

0 pOO 
+ 1 / e 


c(a,a)Fi(a) + c*(a,g.)F2(a) 


da 


-o0 


0* (a, q.) Fa (a) da — e~^ B(a,q.)Fi(a) 


da 


1 Ha r ■ 

/ c(a^q.)Fi(a) + &*(a,q.)Fa(a) 

Q. . / 1 — a I 

j^c*(a,q.)Fa(a) _s(a,q.)Fi(a) 


( 1 — a)0 


q.fi / 1 - a 

o 


da 


da 


(76) 


Further Investigation is required to determine precisely under what 
conditions the flow is 1 aroimd a closed body and uniform, at infinity. 
These investigations seem within areach by the method here indicated. 

This could yield mixed subsonic-supersonic flows with subsonic free- 
stream velocities under the assuiig)tion of the adiabatic eg.uation of state. 

As an application of the method outlined in this paper, fluids of 
which the flows are everywhere subsonic (more precisely, flows satisfy- 
ing the linearized eq.uation of state) are chosen. Precise conditions 
are obtained for flows around closed bodies with prescribed velocities 
at infinity. Flows under these hypotheses have been investigated with 
much success by Chaplygin, von Karm^, Tsien, Bars, and others. Of par- 
ticular Interest in this connection for flows around closed bodies is 
the work of Bers (reference 8) . _ _ 

In the seq.uel to the present report (reference 13) a detailed inves- 
tigation of the flow around a circle is m^e (under the assung>tion of the 
linearized eq.uation of state) . This is done primarily because the flow 
around a circle appears to be basic for the study of flows eiround given 
bodies, as in the case of incompressible fluids. 


CCMIBESSIBLE FLUIE6 UNBEE THE ASSTMFTIOH OF THE LIHEABIZED EQUATION OF STATE 


Of the four basic relations of the flow, eq.uations (l) to (4), the 
first, the equation of state, is replaced by 



( 77 ) 
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where Pi Xs the density for a given velocity qi, and is the cor- 

responding p3?eesure. 

If pressure— volume curves are drawn of equations (l) and (77)^ it 
is ohser'ved that when k is suitably chosen equation (77) represents a 
line tangent ^to tie curve represented by eqmtion (l) at the point pj., 
l/px. Von Knrman and-Tsien used this fact partly to Justify the use of 
the linearized equation of state for the study of subsonic flows. 

From equations (5) euai (77) 


dp p2 


Hence 


a® p® = — k 


From equation (6) 


2 2 

3 ^ + / a_£_ 


if 

2 


dp » constant 


and from equation (79) 


Thus, 


%- + a® p® / ^ = constant 


2 2 

q — a *= consteint 


It follows .that 


q® — a^ 


qS/aS - 1 

l/a2 

M - 1 

-p^/k 

constant 


(78) 

(79) 


(80) 


or 


1 - M" 


= constant 


(81) 
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By ttie normalization, ao * Pq *= 1 > the constants in the right- 

hand sides of equations (80) and (8I) are readily detenained to he minus 
unity eind unity, respectively. Thus, 

a^ =* 1 + -- (82) 


6ind 


From equation ( 77 ) 


p2 = 1 - m2 


» 1 - 


1 + q" 


1 + <l" 


p — Pi =s k 


C/i 


1 __ 


( 83 ) 

( 84 ) 


Again, from the definitions of stream function cp and potential 
function \i/ of an incompresslhle fluid 


also 


dcp = u ds + V dy 
d\^ = — V dx + u dy 


d9 = dx + cpy dy 
d|/ = \1/^ dx + \l/y dy 


Comparing systems (85) and (86), it follows that 


^x “'i'y 
'Py “ -i^x 


(85) 


(86) 


(87) 


These are the Cauchy— Elemann equations. This- estahllshes the well-known 
fact that the cou^ilex potential f2*9+i\l/ of an incon^ressihle flow is 
an analytic fxuiction of the complex variable z = x + ly, the complex 
coordinates in the physical plane. 
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Wlaen the second of e^mtlons (85) is multiplied hy 1 and the two 
equations added 


df5 = (u — Iv) (dbc + i dy) 

-iS 

= qe dz 


(88) 


or 


dz = i an (89) 

q. 

By proceeding as was done in the derivation of equations ( 2 ^), the 
change of the variables (x,y) to the independent variables ( 0 ,q) in 
equations (87) leads to equations 


4 tq 



> 


The transfonnatlon that will symmetrize system 


% = T i (g.) 

1 a (q.) 'i'0 


is given by 



System (90) is therefore symmetrized by the transformation 

1 



and reduces to the Cauchy— Elemann equations 


(90) 


(91) 


(92) 


( 93 ) 




<p^ = -^l/0 


( 94 ) 
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Thus^ the complex potential of an.-inccurpresslhle flow, given in terms of 
the independent variables of the 0, 5~plane is an analytic function of 
the conplex variable w = 0.+ Iq.. - - 


Equations 




3* \i/ 

p 




1 - M-= 
P4 




(95) 


can similarly be symmetrized by the transformation 


<L 



03 


Equations (95) then become 


Ji- y? 

90 = ^ \j/q 


Jl. - 

9;r = 


It 


% 


(97) 


However, under the assuaqition of the linearized equation of state 

y 1 - ^ ^ 

P 


so that equations (97) become the Cauchy— Eiemann equations 


9 e = 1J/5 
'Pq = - t© 


( 98 ) 
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Again, the coiiplex potential of a compreesihle flow under the asBuiuption 
of the linearized equation of state is an analytic function of the cora- 
plex variable w = 0 + iq. The w~plane shall be referred to as the 
distorted hodograph plane. 

Any analytic function of w may be regarded as the con 5 >lex poten- 
tial of either ein incompressible flow in the physical plaine or in the 
hodograph plane, or a con^iressible flow in the distorted hodograph under 
the assumption of this section, and vice versa; then, in general, to 
every incompressible flow around a given body with a given free-stream 
velocity there corresponds a conpressible flow (everywhere subsonic) 
around the same body with the same free— stream velocity. 

This correspondence can be expressed as follows: Given the coaiplei 

potential (w) of an incompressible flow, or of a compressible flow 
in the distoxdied hodograph, then there exists an analytic function 
w = e(0 such that SB is the complex potential in the distorted 

hodograph a given con 5 >reseible flow. 


Since 


and 


where 


Hence 


Jl - if 




and p 




1 + yi + q^ 


( 99 ) 


( 100 ) 


( 101 ) 


gSL- 

1 + yi + qS- 


( 102 ) 
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and 


Also 


i = ^ _ ®_ 

Q 2 ^ 2K 


JL = 

0. 


■4 + - 

2e^ 2 


(103) 


(104) 


Thus, eq.uatlon (I 6 ) may "be put into the more convenient form, 


dz = 


QL 

.10 




(i 


T 10 

1 \ ^ 6 
+ } dT2 + 

2 


’• (- - -) 


dO 


2 ''q. pq.'^ 2 'qt pq.' 

From eq.uations (IO 3 ) and (104), eq.uati6n (IO 5 ) hecomes 


( 105 ) 


e — K e*^ — 

dz = ^ dQ - ^ ® 


3^®. - 

T 


2 K 


dn 


K iw 
= — e 
2 


d:^ 



6n 


(106) 


where w = 6 + iq. Thus, 



where n(w) , the complex potential of a con 5 )ressl'ble fluid in the dis- 
torted hodograph plane, is an analytic function of w. 

Since the con^ilex potential of a ccmgjressl'ble fluid remains the 
complex potential of a coa^ressihle fluid under an analytic transforma- 
tion, set 


v(£) = -1 los 


(108) 
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or 

elw = SflOi (109) 

and 

e-lw . mlltl (110) 

where f(t) and G(t) ^ n[w(t)] are analytic, f($) being arbitrary 
and G( ^) the con^ilei potential of an InconipreBBible flow around a 
given body in the physical ^ -plane. 

When w is considered a function of equation (107) i^iay be 
written as 

z =. {i)6i - ^ (111-) 

When equation (108) is substituted into equation (ill), equation (ill) 
becomes 


z = f«) (112) 

It should be emphasized that f (^) is an arbitrary analytic function of 
while G(^) is an Incongjreasible flew in the f-plane, ^ being 
the physical coordinates' of the flow, and G(^) = r2[w(^)] is the corre- 
sponding cong>ressible flow in the w— plane. The relation. 


Ira G($) = Ira n&f(0] 
= constant 


represents streamlines of the Incouqjressible flow as well -as stre amlin es 
of the compressible flow. If ^ moves along a streamline of the incom- 
pressible flow, then w, • given by equation (108), traverses a stream-- 
llne of the compressible flow in the hodograph plane and z, by the 
treinsformation (108), traces out a streamline in the physical pleine of 
the compressible flow. 

It is convenient to consider G(n the incorapressible flow Eiround 
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a circle of radius R, with circulation and having a free-stream 

velocity namely 

0(t) ^ I (“-3) 

Since e^^ = it follows from equation (llO) that 


Thus, 


. K Glia 

2 f«(^^) 



C ^*(0 

f’(t) 


(114) 


and recalling the relation (102), 


a 

2K e 

K — e ^ 


21 

f*‘ 


1 


1 G^ 
k f « 


(115) 


If, then, f(t) is so determined that a prescribed flow around a 
given body is obtained, it follows that for a given value ^ equation 
(U2) deteimj.nes a point of the flow and equation (115) determines the 
velocity at that point. 

At a stagnation point q = 0. From equation ^115) , q = 0 when 
G'(t) = 0* Thus, stagnation points of the flow occur wherever stagna- 
tion points occur in the inconpressible flow G( ^) in the plane . 

If as I — > oc^ z — > and G’ (^ ) is bounded away from zero at 
Infinity, then, from equation (U5), -f*(0 ™ist be regular and unequal 
to zero at infinity, if the flow is to have a velocity q^ 7 ^ 0 at in- 
finity. Thus, the most general form that f(^) cein have is 

f (U = + ^oU + log t + V l3n ^ 0 (116) 

n=2 
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and 


»(t) 


1 ^ 1 
"bo + "bi > (n - 1)1511 - n 


(117) 


whare the values of h may be complex constants. 


It is known from incongsressible-flow theory that uniform flows past 
closed bodies must also have this form. Therefore, the second term in 
the right-'hand side of eq,uation (112) taJces the form 


PlSlUXtii 

J f'(t) 

“ C_a + Cq I + Cl log 




n -1 


n*e 


(118) 


If the flow at infinity is to be horizontal and of magnitude 
then from eq.uation (ll4) and the fact that 0 


llm I £L(il . 1 

t — 5=- “2 f* ( t) 

When G(z) is given by equation (113), 


(119) 


lim 


f‘(^) 


CZ 



K 

* 2 


(120) 


1 ^ m 

where ^ (t) " ^o Is the condition that the flow at infinity be 

horizontal and q = q , 

03 

In order to ermine the eliape of the body in the stream of the flow, 
the circle t * Ee^'^ is mapped through equation (112) into the z— plane-. 
In order for the circle to map into a closed body, it follows from equa- 
tions (ll6) and (ll8) that 
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Tdi icp ~ Cl icp = 0 


or 


"bi - Cl = 0 


( 121 ) 


It is of Interest to note that, vhen f(0 ~ t, t formula of 

Tsien (reference 6) is ohtalnedj namely 


z 





(122) 


The coefficient in the integral tem ±s different only hecause of a dif- 
ferent normalization. 


Again, when 



Bers’ formula is obtained; namely 


z 




it) 


i+i/n 




(123) 


Here, too, the formulas differ only by a noimalizing factor. 

In Bers* formula n is arbitrary within limits. This freedom 
enables him to determine the conditions for the flow around a closed 
body. Since f('^ in formula (112) is arbitrary and analytic, there is 
an infinite number of arbitrary coefficients. The first two coefficients 
determine the flow at infinity and that the flow be around a closed body. 
The other arbitrary coefficients can be fixed to determine the flew 
around a given body. 

A similar techniq,ue to that developed by Theodorsen and Garriclc 
(references l4 and I5) for determining the coefficients of f(^) might 
be developed in order to obtain a prescribed body. Bers has employed 
this approach by another method with some success. It should be noted 
that when M = 0 the transformation (112) reduces to the initial trans- 
formation used by Theodorsen (reference l4) . This iH5>lies that if an 
integral eq.uation were set up from equation (112) it would reduce to 
Theodorsen’ s when M = 0, so that the integral equation would be a 
generalization of Theodorsen’ s transformation. 
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In another paper Bartnoff determined f(^) such that the circle 
^ goes into a unit circle In the z-plane (accuracy to within 

a few percent), thus giving the oomproesihle flow around a circle. This 
has heen done hy others, notably by Bers. 

The right-hand side of the transfojmiation (112) is invariant uncJer 
an analytic transformation. For, let - 

t •=■ ^(|) ( 12 *^) 

be an analytic function of then by direct subctltutlon ff oq.uation 

(l2h) into eq.uation (112), the trrniafomeitlon becomes 

z = f (125) 

- J J d/d|f[t(|)] 


Because the transformation (112) is invariant under a conformal 
transformation the subsonic flow of a corctproasible fluid around a circle 
is of basic Importance, 


Syracuse University, 

Syracuse, New York, July 1, 19k^. 
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